Hydrodynamic instability during non-uniform growth of a helium crystal 
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We analyze an analog of the hydrodynamic Ray leigh- Taylor instability for the liquid-solid phase 
interface under non-uniform growth of the solid phase. The development of the instability starts 
on conditions of an accelerated interface growth and if the magnitude of acceleration exceeds some 
critical value. The plane and spherical shapes of the interface are considered. The observation of 
the instability can be expected for helium crystals in the course of their abnormal fast growth. 
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Introduction. The fast growth kinetics of the crystal- 
superfluid helium interface results in appearing the insta- 
bilities similar to those inherent in the interfaces between 
fluid media. For instance, it is shown that the crystal sur- 
face is unstable in the steady tangential flow of a fluid 
[l], 0] , demonstrating thus the Kelvin-Helmholtz type of 
instability. The phenomenon is observed qualitatively as 
a distortion of the crystal surface in the fluid jet [§]. The 
onset of such instability occurs when the harmonics in- 
creasing exponentially in time appear in the spectrum of 
crystallization waves. 

It is known after Q that the wave spectrum does not 
vary in the case of the normal incidence of the fluid onto 
the surface. Hence the growth of a crystal should not 
lead to an instability of its surface. However, in the ex- 
periment on the free growth of the crystal initiated at the 
needle point immersed into the liquid bulk 0, [B| one has 
observed an instability at the initial stage of the abnor- 
mal fast growth. After nucleation the crystal grows first 
with the clear kinetic prism-shaped faceting. Then, if 
the initial overpressure at which the crystal is nucleated 
is sufficiently large, the crystal after 0.1-0.2 ms takes the 
shape close to the spherical one with a sizeable ripple on 
the surface, see Fig. [TJ The magnitude of the thresh- 
old overpressure lies within 6-8 mbar and the threshold 
of the instability does not depend upon the temperature 
within errors [6j. 

The attempts to explain this phenomenon within the 
known theoretical works have been unfruitful [6]. The 
essential point, which differs the experiment from the 




FIG. 1: The growth of a crystal at 0.47 K and initial overpres- 
sure 5.2 mbar. The left frame corresponds to 0.19 ms after 
the crystal nucleation. The right frame is taken at 80 ms after 
the crystal nucleation at the pressure close to the phase equi- 
librium pressure. The vertical size of the frames is 2.4 mm. 



cases considered theoretically, is unsteady growth condi- 
tions. Under such large overpressures the crystal grows 
in the oscillating regime with the typical oscillation pe- 
riod ~0.4 ms and reaches the size ~1 mm for the first 
half-cycle. The linear velocity for the motion of the crys- 
tal boundary runs up to ~10 m/s. The interface-directed 
velocity of the liquid phase is one order of the magnitude 
as smaller due to difference in the densities of the liquid 
and solid phases. During the first half-cycle the crys- 
tal growth rate varies from maximum to zero one. This 
gives an estimate ~10 5 m/s 2 for the acceleration of the 
interface and ~10 4 m/s 2 for the fluid beside the crystal. 
Another essential point here is a direction of acceleration 
which is opposite to the velocity of the crystal growth. In 
the non-inertial reference frame related to the interface 
segment this results in the acceleration directed from the 
crystal to the liquid. Since the crystal density is larger 
than the liquid one, the situation proves to be looking like 
the normal Rayleigh- Taylor instability @, Q for the in- 
terface between two classical immiscible fluids of various 
densities. The interface is unstable or stable according 
to whether the acceleration is directed from the lighter 
to the heavier fluid or vice versa. 

However, the fast growth kinetics of a helium crystal 
in the abnormal state gives rise to a series of qualita- 
tive features. The classical Rayleigh- Taylor instability 
develops at the fixed volume of each phase. In contrast, 
the impetuous growth of the liquid-crystal interface can 
result in varying the interface profile due to phase conver- 
sion. In addition, besides the kinetic liquid and interface 
surface energies it is necessary to involve the energy of 
phase conversion depending on the imbalance from the 
phase equilibrium. Thus, in spite of superficial resem- 
blance such instability has specific features which differ 
it from the familiar Rayleigh- Taylor instability. 

The variation of the surface shape of a crystal at its 
non-uniform growth has been studied only in the exper- 
iments on abnormal crystal growth The classical in- 
stability of the crystal-liquid helium interface in the field 
of g ravity is observed by Demaria, Lewellen, and Dahm 
[llj . In these experiments a cell in which the solid and 
liquid phases occupy initially the lower and upper halves, 
respectively, is inverted mechanically by 180°. After in- 
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version the crystal starts to melt, the interface descend- 
ing along the walls and ascending in the centre of a cell. 
The authors gave only the schematic description of the 
process without concrete spatial scale and time refer- 
ence points over technical failures with determining the 
interface position from the filming of the process. Un- 
fortunately, the temperatures at which the experiments 
are performed are not indicated with the exception of 
the crystal growth temperature 1.1 K and maximum one 
1.7 K. The question about the effect of roughening tran- 
sitions upon stability of the crystal facet in the field of 
gravity has remained open. We note only the following. 
The stability observed for the shape of the crystal grown 
at the needle point [4j with its lower facet under condition 
favorable for developing the Rayleigh- Taylor instability 
evidences that the acceleration of gravity did not result 
in any instability of the atomically smooth interface for 
the sizes not exceeding 10 mm at least. 

In the present work we give a theory of the phase in- 
terface instability for the case of the non-uniform inter- 
face growth rate. A comparison is presented between the 
calculation results and the experimental data on devel- 
oping the instability of the crystal shape at the abnormal 
growth. 

The growth kinetics of the crystal surface. Let 

the plane liquid-solid interface grow at the rate V = L(t) 
and regions z > L(t) and z < L(t) be occupied with 
the superfluid and the solid 4 He, respectively. We con- 
sider stability of the growing interface with respect to 
its small distortions from the plane shape, i.e., interface 
z-coordinate is taken as Z = L(t) + £(x,y, t). Assum- 
ing the crystal growth rates to be small compared with 
the first and second sound velocities, we consider the 
hydrodynamics of a superfluid in the approximation of 
incompressible liquid and constancy of the entropy den- 
sity. In this case 0] the equations for the normal and 
superfluid motions can be separated and the potentials 
of the normal and superfluid velocities satisfy, respec- 
tively, V 2 </>„ = and V 2 s = 0. Thus, the solutions 
with wave vector q = [q Xl q y ) parallel to the surface can 
be represented as 

4> s = u s (t)z + A a (t) exp(iqr - qz) 
4>n = u n (t)z + A n (t) exp(iqr - qz) , 

where r = (x, y) and velocities u s and u n correspond to 
those of the undisturbed motion. The pressure in the 
fluid is a sum of pressures P — P n + P s and 

P S =P SOO - p s {<j>s + (V</g 2 /2) + p.giLoo - z) 

Pn = Pice - Pn(4>n + (V0„) 2 /2) + /9„g(ioo ~ z) . 

Let us employ the boundary conditions used in the 
theory of crystallization waves and crystal growth kinet- 
ics. From the continuity of the mass flow we have at the 
interface z = Z(x, y, t) 

jy = PnV n v + p s V s „ = (p - p')Z . 



Here v is the normal to the surface. 

To simplify the consideration, farther on we neglect 
shear components a^k for the stress tensor <7j& in the 
solid phase. In other words, we suppose the stress tensor 
to be diagonal, i.e., c^. = —P'S^. Then the solid looks 
like a liquid under pressure P 1 . 

From the continuity of the momentum flux density and 
involving the surface tension and interfacial curvature, we 
arrive at 



P 



Zuf + Ps {v a - Zu) 2 
= ^ikdZ 2 /dxidxk , 



(P' + p'Z 2 ) 



where 7^ = otSik + da 2 / d(fiidipk is the surface stiffness 
tensor expressed in terms of surface tension a and deriva- 
tives over angle if between the normal to the displaced 
surface and the axis z. The quantities with dash are 
referred to the solid phase. 

For the normal velocity according to work [lfj| , we sup- 
pose the sticking of the normal component to the surface 
as for a viscid fluid 

Vnv = Z . 

The growth rate of the solid phase is described by the 
kinetic growth coefficient K and connected with the dif- 
ference in chemical potentials per unit mass as usual 



Z = K 



(v s - Zu) 



P + T 



Then we have for the difference in the chemical potentials 
and temperature Q 
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where a is the entropy and the quantities with index "00" 
stand for the magnitudes taken far from the interface. 
Pressure P c denotes the equilibrium pressure at which 
the phase transition takes place. 

Solving a set of the boundary conditions, one finds 
that the undisturbed motion of the plane interface at 
rate V = Lit) obeys 



V— 

K 



P -P 



[AP + pg(L 00 -L)]+pJVL+ — 
Here p e f is the effective density of the interface 

Pef = Pn + (/?' - Ps) 2 /Ps 

and AP = P^ — P c is the overpressure. The equation for 
the disturbed motion in the linear approximation reads 

Pel - + 77 C + [nkqiqk + (p - p)g - p C fV] ( = 0. 

q K 
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For the uniform growth of the interface V = 0, the equa- 
tion goes over into the familiar dispersion relation for the 
crystallization wave spectrum 13, 

The instability. Let us consider the stability of the 
flat-shaped interface with respect to small perturbations 
£ ~ exp(At) for the uniformly accelerated growth of the 
crystal. The root with ReA(g) > corresponds to insta- 
bility. For the acceleration exceeding the threshold one 
V c = g(p' — p)/ p c i, the interfacial distortion will increase 
for the wave vectors satisfying jikqiqk < PdV — (p' — p)g. 
Note that V c does not depend on the growth coefficient, 
i.e., on the dissipative properties of the interface, and is 
positive. The latter corresponds to the case when the 
interfacial acceleration is directed to the fluid. Thus, the 
stability due to non-uniform growth of the plane interface 
appears only for the accelerated growth of a crystal. 

The value qo corresponding to the maximum magni- 
tude Re X(q) meets the shortest time for the development 
of the instability which will be characterized by the spa- 
tial scale 27r/<2o- The value qo can be found from the 
equation 



2 = q\ _ V2 J_ 
% 3 3 JjpVt K 



3/2 



where q c is the value related to the upper bounds of in- 
stability according to jq 2 = p ci V — g{p' — p). For the 
large magnitudes of the growth coefficient or large accel- 
eration while K 4 V ^> p lA /"fp% 1 values qo and \(qo) are 
equal approximately to 



qo 



An 



gc 

V3 



= ( ^ q o 

\3V3 Pa 



(6V3~fp cl K 2 q c y/ 2 

V2 „ , 



2\/3 PdK ' 

In the opposite limit K 4 V <C p' 4 /7/°cf one has roughly 



1/3 



, 7Pcf r , 2 4 

A = j jq c 



On the whole, values qo and Ao decrease as the ki- 
netic growth coefficient reduces or dissipation of the in- 
terface enhances. From the experimental point of view 
this may require the crystal surface of sufficiently large 
sizes d > 2w/qo and the large time of supporting the 
accelerated growth t > 1/ Ao in order to realize the inter- 
facial instability at the uniformly accelerated growth. 

The spherical geometry. In the experiment a crys- 
tal nucleates at the needle point and then grows free. The 
anisotropy of the kinetic growth coefficient is not large. A 
ratio of the maximum to minimum growth rate does not 
exceed 2-3. Let us consider the stability of the spherical 
shape of the growing crystal. For simplicity, we assume 



the isotropy of the surface tension a. As it follows from 
the estimates given in the Introduction, the acceleration 
of gravity is negligibly small as compared with the typical 
magnitudes of the interfacial acceleration. 

The equation of the interface is given as r = R s (t, f2) = 
R(t) + C(t, fi) where C = Q(t)Yi(Q,) is a fluctuation of the 
interface expanded in the spherical harmonics of degree 
I — 0, 1, 2 . . .. We seek for the velocity potentials of the 
normal v n = V</>„ and superfluid v s — V0 S motions in 
the form 

4> s = ~u s R 2 /r +A l Y l /r l+1 
K = -u n R 2 /r + BiYi/r^ 1 . 

Using the same boundary conditions as above, we find 
the velocities of the undisturbed flow of the liquid phase 

u s = R (p s - p')/p s , u n = R, 

and coefficients Ai and Bi describing the disturbed mo- 
tion of the interface 



A, = - 
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Employing relation for the pressure and the dependence 
between growth rate and chemical potential difference, 
we obtain for the undisturbed growth of the solid phase 



iK.iRR+^R K 



2a 



p'-p R 



The growth equation looks like the motion of a parti- 
cle with the effective mass M(R) — 47rp fi? 3 , drag force 
4ttR 2 p'K~ 1 R, and potential energy U(R) = 4iraR 2 - 
((p'-p)/ P )AP(4^ 3 /3). 

The small deformations for the spherical surface of the 
solid phase are described by the relation 



l + l 



l + l 



K 



a ^ p e {——R)Q=0. (1) 



R 2 



l + l 



For R = R = and in the lack of energy dissipation, the 
equation describes the spectrum of crystallization waves 
on the spherical interface 



PcfR 3 



(l-l)(l + !)(/ + 2) 



The asymptotic behavior of function Q (t) at large t is 
given by 



Ci(*)~exp dt' 



-W) + ^Jx 2 (t') + x t (t')-^(t') 



FIG. 2: The phase diagram for the stability of a spherical 
crystal at the non-uniform growth. The dashed parabola sep- 
arates the regions of oscillating (below the curve) and relax- 
ation modes. The regions of instability are shown with shad- 
ing. The curved arrow demonstrates schematically, but in 
the scale, the trajectory of the crystal growth during the first 
half-cycle of the oscillating growth in the abnormal state. 

where 

Xi(t) = h - — , 

w 2 R 2 p c f KR 

ui(t) = lj 2 01 -(1-1)R/R. 

For the fixed time instant, the proper frequencies are 
determined by equation (JTJ . The onset of instability at 
the fixed time occurs if the imaginary part of the fre- 
quency becomes positive. To estimate, we consider the 
perturbation as Q ~ exp(— iwi). In the dimensionless 
units 

w , 1 / 3R p' l + l-l R 

Ct= , £ = 1- — i V = — o 

tool 2w ; V R p c f KR I R 

the proper frequencies are determined from 

The condition of appearing the instability is Im(f2) > 0. 
The result is given in Fig. [2] Provided that 1 — r\ — £ 2 > 0, 
the frequencies have the real part but the instability ap- 
pears for the negative value of parameter £. If the radi- 
cand is negative, the both frequencies are imaginary. In 
this case the left half-plane of the phase diagram corre- 
sponds also to instability. The appearance of the insta- 
bility takes place in the right half-plane for r\ > 1. 

In Fig. [3] a simulation is shown for the behavior of 
spherical harmonics at the initial stage of the abnormal 



FIG. 3: The damping of the spherical harmonics during crys- 
tal growth as a function of the time and degree of harmonic. 
The numbers beside the lines denote the degree of harmonic. 
The normalized time t = 2nfoT is put along the abscissa axis. 



crystal growth. The variation of the crystal radius is de- 
scribed approximately by expression R = Rq sin 27t/ot. 
The following values Rq = 0.1 cm, K = 0.1 s/cm, 
fa = 2.5 kHz of experimental parameters at 0.5 K and 
overpressure 6 mbar are taken to calculate. The other 
quantities are the tabular data. The deformation of unit 
amplitude is put for the initial values. The velocity is 
varied within (-100, 100). In the figure the solutions are 
given with zero initial velocity. Note that the variation 
of the initial velocity in the range analyzed does not lead 
to the qualitative change in the behavior of the curves 
with the exception of first move at which fast relaxation 
occurs. Next in all cases the damping of the harmonics, 
checked to I = 200, is observed. The increase of the de- 
gree of the harmonic, as it follows from the simulation, 
results in the faster relaxation of the deformation. Thus, 
the numerical simulation agrees with the previous con- 
clusion about the crystal interface stability at the initial 
growth stage. 

Summary. Non-uniformity of the crystal growth re- 
sults in new qualitative hydrodynamic instability of the 
liquid-solid interface. At the plane interface the develop- 
ment of instability starts only for the accelerated growth 
of the solid phase when the magnitude of the interfacial 
acceleration exceeds the critical one determined by the 
acceleration of gravity. 

For free growth of the solid phase when its shape is 
almost spherical, the conditions of appearing the insta- 
bility are determined both by the interfacial acceleration 
and by the growth rate. The experimental situation in 
which one may expect an observation of such instability is 
realized with the abnormal helium crystal growth when 
the interfacial mobility increases drastically by several 
orders of the magnitude [l3| and the interfacial velocity 



and acceleration reach the huge magnitudes as is noted 
in the Introduction. For lowest mode 1 = 2, the normal- 
ized order-of-magnitude parameters are about £ ~ 10 2 
and r\ ~ — 10 4 . However, during the first half-cycle of the 
oscillating growth, as is seen from Fig. [51 the trajectory 
of the crystal growth lies in the fourth quadrant in which 
the instability considered here does not appear. Thus the 
involvement of the non-uniform interfacial motion alone 
may fail to explain the liquid-solid interface instability 
observed. Note that none of the hydrodynamic instabil- 
ities known so far can be responsible for the instability 
of the crystal shape observed at the abnormal crystal 
growth [6J]. 

Acknowledgments. The work is supported by the 
RFBR Grant No. 05-02-16806. 



* burmi@kurm.polyn.kiae.su 



5 



[1] M. Yu. Kagan, Zh. Eksp. Teor. Fiz. 90, 498 (1986) [Sov. 

Phys. JETP 63, 288 (1986)]. 
[2] P. Nozieres, M. Uwaha, J. de Physique 47, 263 (1986). 
[3] L. A. Maksimov, V. L. Tsymbalenko, Zh. Eksp. Teor. 

Fiz. 122, 530 (2002) [JETP 95, 455 (2002)]. 
[4] V. L. Tsymbalenko, Fiz. Nizk. Temp. 21, 162 (1995) [Low 

Temp. Phys. 21, 120 (1995)]. 
[5] V. L. Tsymbalenko, J. Low Temp. Phys. 121, 53 (2000). 
[6] V. L. Tsymbalenko, Zh. Eksp. Teor. Fiz. 126, 1391 (2004) 

[JETP 99, 1214 (2004)]. 
[7] G. I. Taylor, Proc. Roy. Soc. (London) A201, 192 (1950). 
[8] G. Birkhoff, Hydrodynamics, (Princeton Press, 1960). 
[9] L. D. Landau, E. M. Lifshits, Fluid Mechanics, (Perga- 

mon Press, 1987), §140. 
[10] A. F. Andreev, V. G. Knizhnik, Zh. Eksp. Teor. Fiz. 83, 

416 (1982) [Sov. Phys. JETP 56, 226 (1982)]. 
[11] C. D. Demaria, J. W. Lewellen, A. J. Dahm, J. Low 

Temp. Phys. 89, 385 (1992). 
[12] S. Balibar, H. Alles, A. Ya. Parshin, Rev. Mod. Phys. 77, 

317 (2005). 

[13] V. L. Tsymbalenko, J. Low Temp. Phys. 138, 795 (2005). 



